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We report Monte Carlo simulations of critical dynamics far from equilibrium on a perfect and 
non-perfect surface in the 3d Ising model. For an ordered initial state, the dynamic relaxation of 
the surface magnetization, the line magnetization of the defect line, and the corresponding suscep- 
tibilities and appropriate cumulant is carefully examined at the ordinary, special and surface phase 
transitions. The universal dynamic scaling behavior including a dynamic crossover scaling form 
is identified. The exponent /3i of the surface magnetization and (32 of the line magnetization are 
extracted. The impact of the defect line on the surface universality classes is investigated. 



PACS numbers: 64.60.Ht, 68.35.Rh, 05.20.-y 



I. INTRODUCTION 

The breakdown of space and time translation invari- 
ance leads to geometric and temporal surface effects. The 
former is very common in a system whose spatial correla- 
tion length is comparable to its dimensions. Such effects 
become even more important when nano-scale materials 
are concerned. In a recent experiment, for example, an 
anomalous temperature profile of the phase transitions 
was observed in the presence of a ferromagnetic surface 
The latter occurs in a nonequilibrium system, which 
is prepared by suddenly quenching the system to its crit- 
ical temperature from any given initial condition. 

The breakdown of space translation invariance modi- 
fies the critical behaviors near geometric surface and new 
critical exponents must be introduced. There may ex- 
ist several universality classes in one bulk system in the 
presence of free geometric surface. The critical behavior 
of geometric surface has been extensively studied, and 
the equilibrium phase diagram has been well established 
in the past decades 0, H, [j, [1] . However, most previous 
studies concentrated on the static behaviors [ 5 [3 , M [ p . fiot 
and the dynamics in the long-time regime [111 1 12l. Il3l| . i.e. 
system only with geometric surface. The critical dynam- 
ics of surface in the macroscopic short-time regime, i.e., 
when the system is still far from equilibrium, is much less 
touched [1]. 

On the other hand, for a system quenched to its crit- 
ical temperature, because there is no characteristic time 
scale, the temporal surface has long-lasting effect. This 
effect has very important consequences. One is that in 
nonequilibrium dynamic relaxation of magnetization, if 
in the initial state, there is small, nonvanishing magne- 
tization Too ^ 1, the magnetization grows as TOpt^ with 
9 being a new nonequilibrium dynamic exponent [T^ . In 
such short-time critical dynamics, there exist two com- 
peting nonequilibrium dynamic processes. One is the do- 
main growth with scaling dimension Xi and the other is 
the critical thermal fluctuation with scaling dimension 



X — j3/v. Because the spatial correlation length ^ grows 
as i^/^, we can relate 9 to Xi and a; by 6* = . Gener- 
ally Xi is larger than x and the net effect is the domain 
growth in the nonequilibrium relaxation process. This 
short-time critical dynamics of bulk has been established 
in the past decade, and successfully applied to different 
physical systems [ll, [H, [H, [13, EE iSl • Based on the 
short-time dynamic scaling, new techniques for the mea- 
surements of both dynamic and static critical exponents 
as well as the critical temperature have been developed 
[20l . [2ll.l23l. Recent pro gresses can be found partially in 
Refs^lM H MMi ■ 

Obviously, the physical phenomena are more compli- 
cated, when both temporal and geometric surfaces are 
considered. The interplay between both surfaces em- 
braces many interesting physics and is worth for care- 
ful studiesjl, 113] ■ Recently it is reported that in non- 
equilibrium states, the surface cluster dissolution may 
take place instead of the domain growth [H, [2^. In 
these studies, the dynamic relaxation starting from a 
high-temperature state is concerned. 

The impact of defect on geometric surface is also of 
great concern. The presence of imperfection may alter 
the surface university classes and even the phase diagram. 
The former is easily signaled from the non-equilibrium 
dynamics as in the case of bulk [2^ [25l[30|. 

In this paper, we study the short-time critical dynam- 
ics on a perfect and non-perfect surface with Monte Carlo 
simulations. We generalize the universal dynamic scaling 
behavior to the dynamic relaxation at geometric surfaces, 
starting from the ordered state. At the ordinary, spe- 
cial and surface phase transitions, the dynamic scaling 
behavior of the surface magnetization, susceptibility and 
appropriate cumulant are identified. The static exponent 
[3i of the surface magnetization and [32 of the line magne- 
tization of the defect line are extracted from the dynamic 
behavior in the macroscopic short-time regime. The ro- 
bustness of surface university class against extended de- 
fect is investigated by means of non-equilibrium dynam- 
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ics. The surface transition and special transition can also 
be detected from the short-time dynamics. 

The remaining part of this paper is organized as fol- 
lows. In Sec. ini the definition of the model and the 
short-time dynamic scaling analysis are presented. In 
Sec. IIIII and lTVl the dynamic relaxation on a perfect and 
non-perfect surface is studied. In Sec. |Vl the results are 
summarized. 



II. MODEL AND DYNAMIC SCALING 
ANALYSIS 

A. Model 

The Hamiltonian of the 3d Ising model with Glauber 
dynamics and line defect on free surface in the absence of 
external magnetic field can be written as the sum of bulk 
interactions, surface interactions and line interactions, 

Z^<y> '^xyzCx'y'z' , 

(1) 

where spin a can take values ±1 and < xyz > indicates 
the summation over all nearest neighbors. The first sum 
runs over all links including at least one site that does 
not belong to the surface, whereas the second sum runs 
over all surface links excluding the links that both sites 
are inside the defect line. The last summation extends 
over all links which belong to the defect line. Jb, Js and 
Ji are the coupling constants for the bulk, surface and 
defect line respectively. For ferromagnetic materials, Jb 
and Js are positive. It is generally believed that the dy- 
namic universality class of Glauber dynamics is insensi- 
tive to the detailed algorithm used as long as the updat- 
ing algorithm is local. Here we use Metropolis spin-flip 
algorithm. Without explicitly specified, the dynamic ex- 
ponent refers to the Ising model with Glauber dynamics 
in the following discussions. 

For a perfect surface, i.e., J; = Js, it is well known 
that there exists a special threshold r^p = Js/ Jb in equi- 
librium. For Js/ Jb < fsp, the surface undergoes a phase 
transition at the bulk transition temperature Tb, due to 
the divergent correlation length in the bulk. This phase 
transition is called the ordinary transition, and the crit- 
ical behavior is independent of Js/ Jb- See Fig. [TJ This 
is a strong universality. For Js/ Jb > fsp, the surface first 
becomes ferromagnetic at a surface transition tempera- 
ture Ts > Tb, while the bulk remains to be paramagnetic. 
If the temperature is further reduced, the bulk becomes 
also ferromagnetic at Tb- The former phase transition is 
called the surface transition and the latter is called the 
extraordinary transition. It is generally believed that the 
surface transition belongs to the universality class of the 
2d Ising model 0,11. Around r^p occurs the crossover 
behavior. At exactly Js/Jb = ^sp, the lines of surface 
transition, ordinary transition and extraordinary tran- 



sition meet at this multicritical point with new surface 
exponents. The surface and bulk become critical simul- 
taneously at this point and this phase transition is called 
the special transition. The best estimate of rsp for the Zd 
Ising in equilibrium is 1.5004(20) [31j. 

For a non-perfect surface, we introduce a defect line 
with coupling strength Ji onto the surface. Generally 
speaking, the impact of imperfection on a surface is two 
fold. Take a surface with random bond disorder as an 
example. The randomness may reduce the surface transi- 
tion temperature, and alter the global phase diagram of a 
semi-infinite system. For example, the special transition 
point of the Ising model with a amorphous surface is lo- 
cated at Tsp = 1.70(1) [111, noticeably larger than that of 
the Ising model with a perfect surface rsp = 1.5004. An- 
other effect of the randomness is that it may change the 
universality class of the surface. The relevance or irrele- 
vance of random imperfections on the pure surface can be 
assessed by the Harris- type criterion j33| . The extended- 
Harris criterion states that for a surface with random 
bond disorder, the disorder is relevant for an > but 
irrelevant for an < 0. Based on this criterion, the ran- 
dom surface coupling of the Ising model is irrelevant at 
the ordinary transition since an < 0. In this case, it 
was rigorously proved fif"^ — /JJ"''' by Diehl based on 
the Griffiths-Kelly-Sherman inequality, where (3f'^'^ is the 
critical exponent at the ordinary transition on a random 
bond surface [3^. The situation is less clear at the spe- 
cial transition, for an is very close to 0. Recent simu- 
lations suggest that an < and hence the disorder is 
irrelevant [lo] . The irrelevance at the special transition 
has also been reported in Ref. [12] ■ At the surface tran- 
sition, the surface is equivalent to the 2c? Ising model. 
The disorder only leads to logarithm correction (see (30| 
and reference therein). In the case of defect line, the 
defect doesn't shift the transition temperatures of the 
surface transition, and therefore, the special transition 
point rsp at which the surface transition line and ordi- 
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FIG. 1: Schematic phase diagram for the semi-infinite Ising 
model with bulk coupling Jh and surface coupling ratio Ja/ Jb- 
Tb is the bulk transition temperature and the ferromagnetic 
bulk is denoted by FB while the paramagnetic is denoted by 
PB. The surface phases are labeled FS for a ferromagnet and 
PS for a paramagnet. 
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nary transition line meet remains to be unchanged [35|. 
We only consider the robustness of the ordinary, special 
and surface transition in the presence of line defect. 

B. Dynamic scaling analysis 



TABLE I: The bulk critical temperature and critical expo- 
nents of the 3d Ising model. 







Z3d 


4.5115248(6) [36] 


0.6298(5) [37] 


2.042(6) [38] 



For a dynamic system, which is initially in a high- 
temperature state, suddenly quenched to the critical 
temperature, and then released to the dynamic evolu- 
tion of model A, one expects that there exist univer- 
sal scaling behaviors already in the macroscopic short- 
time regime ^14j . This has been shown both theoreti- 
call y a nd numeri cally in a variety of statistical systems 
0,111, [IE 113, Ella US, Hi, and it explains also the spin 
glass dynamics. Furthermore, the short-time dynamic 
scaling behavior has been extended to the dynamic re- 
laxation with an ordered initial state or arbitrary initial 
state, based on numerical simulations [16, 23, 25, 38, 39]. 
Recent renormalization group calculations also support 
the short-time dynamic scaling form for the ordered ini- 
tial state [26j . 

On the other hand, Ritschel and Czcrncr have gener- 
alized the short-time critical dynamics in a homogenous 
system to that in an inhomogeneous one, i.e., the systems 
with a free surface, and derived the scaling behavior of 
the magnetization close to the surface for the dynamic 
relaxation with a high-temperature initial state [27[ . Re- 
cent development can be found in Refs. [28l. [2^. In this 
paper, we alternatively focus on the dynamic relaxation 
with the ordered initial state, and with a non-perfect 
surface. As pointed out in the literatures [H, HSl, the 
fluctuation is less severe in this case. It helps to obtain a 
more accurate estimate of the critical exponents at sur- 
face. From theoretical point of view, it is also interesting 
to study the dynamic relaxation with the ordered or even 
arbitrary initial state. 

Similar to the scaling analysis in bulk [3, [H, [H, [11] , 
we phenomenologically assume that, for dynamic relax- 
ation with ordered initial state, the surface magnetization 
decays by a power law, 

<mi(i) >-t-'^i/'^="% (2) 

after a microscopic time scale tmic- Here < • > represents 
the statistical average, f3i is the static exponent of the 
surface magnetization, is the static exponent of the 
spatial correlation length, and Zs is the dynamic expo- 
nent. This assumption can be understood by noting that, 
for nonequilibrium preparation with ordered initial state 
niQ = 1, the dynamic relaxation is governed by critical 



thermal fluctuation with scaling dimension xi = Pi/vs- 
For the ordinary and special transitions where the crit- 
icality of surface originates from the divergence of the 
correlation length in bulk, there are no genuine new sur- 
face dynamic exponent Zs and static exponent . Vs and 
Zs are just the same as those in the bulk, i.e. Vs = fad 
and Zg = Z3d, while (3i is neither that of the 2d Ising 
model nor that of the 3c? Ising model [ll|. For the sur- 
face transition where the critical fluctuation of surface 
is of the universality class of the 2c? Ising model, it is 
generally believed that all static and dynamic exponents 
are the same as those of the 2d Ising model 0,11]. i.e. 
Ih = f32d = 1/8, i^s = V2d = and z^ = Z2d « 2.16(2) 

Another important observable is the second moment 
of the surface magnetization, or the so-called time- 
dependent surface susceptibility, defined as 

Xii = L^[< mj > - < mi >2]. (3) 

Simple finite-size scaling analysis reveals that 

Xii(t)^F"/''=^^ (4) 

Here the exponent ju/vs is related to (ii/v^ ^n/vs = 
d—1 — 2(3i/i's, with c? = 3 being the spatial dimension of 
bulk. This is nothing but the scaling law in equilibrium 
between the exponent of the surface susceptibility and 
the exponent of the surface magnetization. One can also 
understand the scaling behavior in Eq. @ in an intuitive 
way. In equilibrium, xii behaves as xii ~ LPi^^l"' with 
L being the lattice size. In the dynamic evolution, X\\(t^ 
should be related to the non-equilibrium spatial correla- 
tion length ^(i) with xii(i) ~ since the finite 
size effect is negligible. Then the growth law ^(t) '-^ t^l^" 
of the non-equilibrium spatial correlation length imme- 
diately leads to Eq. (H]). 

Alternatively, one can also construct the appropriate 
time-dependent cumulant C/(t) =< > / < mi >^ —1. 
Obviously, \J{t) - i(7ii+2/3i)/i^.z. ^ Ytqtc the scahng law 
lii/^s = d—^ — 2(3i/i>s, one derives (711-f 2/3i)/i^s = c?— 1. 
The scaling behavior of U{t) then reduces to the standard 
form ^16i. i38>]. 

C/(t) ~ it'^-i)/^^, (5) 

with d—1 being the spatial dimension of the surface. 

In other words, from Eqs. ^ and ([1]), or from Eqs. ^ 
and ([5]) , we obtain independent measurements of two crit- 
ical exponents, e.g., Pi/vs and Zg. Alternatively, if we 
take Vs and z^ as input, we have two independent esti- 
mates of the static exponent (3i of the surface magneti- 
zation. This may testify the consistency of our dynamic 
scaling analysis. 

All foregoing equations involve the bulk exponents Vg 
and Zs ■ Therefor an accurate estimate of the surface crit- 
ical exponents /3i needs precise values of ly^d and 23^, as 
well as Z2d- Since the 3c? bulk Ising model has been ex- 
tensively studied with various methods, many accurate 
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TABLE 11: Summary of the surface critical exponents at the ordinary and special transition in the 3d Ising model, as obtained 
by different techniques. MF: mean-field, MC: Monte Carlo simulations, FT: field-theoretical methods, CI: conformal invariance. 
The data marked with * are calculated by using scaling law 2/3i -I- 711 = (d — f )z/s- 





MF [2] 


MC [8] 


MC [9] 


MC [40] 


MC [fO] 


MC f4f] 


MC-fCf [36] 


FT [42] 


this work 


fjord 
Pi 


f 


0.78(2) 


0.807(4) 


0.80(f) 


0.796(f) 




0.798(5) 


0.796 


0.795(6) 


P? 


1/2 


0.f8(2) 


0.238(2) 




0.229(f) 


0.237(5) 




0.263 


0.220(3) 


ill 


1/2 


0.96(9) 


0.788(f) 




0.802(3)* 


0.785(ff)* 




0.734 


0.823(4) 



results of the critical exponents and transition tempera- 
ture are available. We concentrate our attention to the 
surface exponents and take the bulk exponents as input. 
The results of the bulk exponents of the 3c? Ising model 
are summarized in Table ID The criteria to choose those 
values are their relative accuracy, as well as the methods 
used to extract these exponents. 



C. Simulations 

In this paper, with Monte Carlo simulations we study 
the dynamic relaxation of the Zd Ising model on a perfect 
and non-perfect surface at the transition temperature, 
quenched from a completely ordered initial state. The 
standard Metropolis algorithm is adopted in the simula- 
tions. In order to investigate the surface critical behav- 
ior, we apply the periodic boundary condition in the xy 
plane and open boundary condition in the z direction to 
the L X L X L cubic lattice. 

The main results are obtained with the lattice size L — 
128 and L = 80, and additional simulations with other 
lattice sizes are also performed to study the finite-size 
effect. For a perfect surface, the surface magnetization is 
defined as 

1 ^ 

mi = ^"^{(J^yl+Cr^yL), (6) 

xy 

and its critical exponent is denoted by For a non- 
perfect surface, the defect line is placed at surface posi- 
tion X = L/2 and the line magnetization is defined as 

1 ^ 

xy 

and its critical exponent is denoted by P2 ■ The spin axyz 
denotes the spin sitting at site {x,y,z). We measure the 
surface and line magnetization during the nonequilibrium 
relaxation. We average from 5000 to 20000 runs with dif- 
ferent random numbers to achieve a good statistics. Er- 
ror bars are estimated by dividing the total samples into 
two subgroups, and by measuring the exponents at dif- 
ferent time intervals. Most of the simulations are carried 
out on the Dawning 4000A supercomputer. The total 
CPU time is about 3 node-year. 



III. SHORT-TIME DYNAMICS ON A PERFECT 
SURFACE 

In this section we study the nonequilibrium critical dy- 
namics on a perfect surface, i.e. Ji = Jg- To investi- 
gate the critical behavior on the surface, it is important 
to know the special transition point Vsp- For a perfect 
surface of the 3d Ising model with ferromagnetic interac- 
tions, there exist rather accurate estimates of r^p in equi- 
librium, e.g., rsp = 1.5004(20) in Ref. [31]. We adopt this 
value as the special transition point. As illustrated later, 
the special transition point also be estimated from 

the scaling plot of a dynamic crossover scaling relation. 

For the ordinary phase transition, the dynamic relax- 
ation of the surface magnetization with different Jg/ Jb 
are shown in Fig. [51 The curves of Js/Jb = 1-0 with 
L = 40 and L = 80 overlap up to t > 300MCS (Monte 
Carlo sweep per site). It confirms that the finite-size 
effect is negligibly small for L = 80 up to at least 
t = lOOOMCS, since the correlating time of a finite sys- 
tem increases by ti, ^ L^. In Fig. [21 a power-law be- 
havior is observed for all Js/Jb- The microscopic time 
scale tmic, after which the short-time universal scaling 
behavior emerges, in other words, after which the cor- 
rection to scaling is negligible, gradually increases as the 
surface coupling is being enhanced. For Js/Jb = 0.2, 
tmic ~ lOMCS, while for Js/Jb = 1.2, tmic ~ lOOMCS. 




10 100 t 1000 



FIG. 2: Dynamic relaxation of the surface magnetization is 
displayed with solid lines on a double-log scale, at the ordi- 
nary transition with various Js/Jb, and at the special tran- 
sition Js/Jb ~ rsp = f.5004. The temperature is set to the 
bulk critical temperature Tc, and the lattice size is L = 80. 
Open circle are the data for Js/Jb = f-O and L — 40. Well 
away from the special transition, the slope of the curves is 
independent of Js/Jb- 
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50 t 500 1000 



FIG. 3: Dynamic relaxation of the surface susceptibility and 
cumulant at the special transition is plotted on a double-log 
scale. The lattice size is L = 128 and T = Tc. 



A direct observation in Fig. [5] is that the curves for 
Js/Jb < Tsp are parallel to each other. By fitting these 
curves to Eq. ©, we obtain /3f = 0.790(7), 0.792(6), 
0.795(6), 0.786(6) and 0.755(12) for J,/Jb = 0.2, 0.5, 
0.8, 1.0 and 1.2 respectively. The values of /JJ"''^ at 
Js/Jb = 0.2, Js/Jb — 0.5 and Js/Jb ~ 0.8 are well con- 
sistent with each other within error. It indicates that 
the ordinary transition is universal over a wide range in 
the Js/Jb space. Deviation occurs for Js/Jb > 1-0 and 
manifests itself as the effect of the crossover to the spe- 
cial transition. This is in agreement with the observa- 
tion in Ref. [8)]. From our analysis, /3f = 0.795(6) is a 
good estimate for the ordinary transition. In Table HIl we 
have compile all the existing results which were obtained 
with simulations and analytical calculations in equilib- 
rium, and our measurements from the non-equilibrium 
dynamic relaxation. A reasonable agreement in /3j"''' can 
be observed. Part of the statistical error in our dynamic 
measurements is from the input of the bulk exponents i^s 
and Zs- 

With nii{t) at hand, one may proceed to investigate 
the time-dependent susceptibility, Xii(^)- the case of 
the ordinary transition of the 3d Ising model, however, 
7ii is negative. Therefore the xii is suppressed during 
the time relaxation according to Eq. ^ and fluctuating 
around if nonequilirium preparation is an ordered ini- 
tial state Xii(O) = 0- The power-law behavior in Eq. (|4]) 
could not be observed. Nevertheless, at the special tran- 
sition, where 711 is positive, the situation is different. 
The power-law behavior of the surface susceptibility and 
cumulant shows up. 

In Fig. [2] and [3l the surface magnetization, surface 
susceptibility and appropriate cumulant are displayed at 
the special transition Js/Jb = ^sp. A power-law behav- 
ior is observed for all three observables. From the slope 
of the curve of the surface magnetization, we measure 
Pl^/vsZs = 0.171(2), and then obtain = 0.220(3) 
with Vs and Zg in Table U as input. From the curve of the 
surface susceptibility, we measure jli/iygZs — 0.640(3), 
and then calculate 7^^ — 0.823(4). From the scaling law 
jii/i^s = d - 1 - 2/3i/z/s, one derives = 0.218(2), 
which is in good agreement with /3j^ = 0.220(3) esti- 
mated from the surface magnetization. The scaling be- 
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FIG. 4: The scaling plot of mi{t) according to Eq. ([9]) around 
the special transition Js/Jb = fsp = 1.5004. The time window 
in this plot is within [10, 1000]. The lattice size is L = 80 and 
T = Tc 

haviors in Eqs. ^ and ^ indeed hold. 

The remarkable feature of the cumulant on the surface 
is that its scaling behavior in Eq. ([5]) does not involve 
the exponent (3i of the surface magnetization. From the 
curve in Fig. [31 we obtain (d — l)/zs = 0.996(11), then 
calculate the few/fc dynamic critical exponent Zs — 2.01(2). 
This value of Zs is very close to z^d — 2.04(1) measured 
in numerical simulations in the bulk in Table [H and it 
confirms that the dynamic exponent on the surface is 
the same as that in the bulk. 

In order to describe the dynamic behavior of the sur- 
face magnetization around Tgp, we need to introduce a 
crossover scaling relation. To understand the scaling rela- 
tion in non-equilibrium states, we first recall the crossover 
scaling relation in equilibrium. In equilibrium, toi(t) 
near the special transition is described by a crossover 
scaling relation 

mi(T)r-^i'' = M,g{T-^{Js/Jb - rsp)), (8) 

where r = 1 — T/Tc is the reduced temperature, and (j) is 
the crossover expoirent. From the crossover scaling rela- 
tion of toi(t) , one can determine the special transition 
point rsp as well as and (j) Nevertheless, up to now 
it has not been studied whether there also exists a cor- 
responding crossover scaling relation in non-equilibrium 
states. Here we will verify that such a dynamic crossover 
scaling form indeed exists. For simplicity, we consider 
the case when Js/Jb approaches the special transition 
from rjp and the system is at the bulk critical tempera- 
ture. Now the non-equilibrium spatial correlation length 
^(t) ~ t^/^ takes the place of the equilibrium spatial cor- 
relation length r^"^. By substituting t"^/'^"^" for t into 
Eq. ([5]), we obtain 

^^(i)i/3r/-.-. = M,,,g{t^/''^^^ {Js/Jb - rsp)). (9) 

We have performed non-equilibrium simulations at 
Js/Jb = 1.30, 1.35, 1.37, 1.40, 1.43, 1.45, 1.47, 1.49, and 
made a scaling plot according to Eq. © . This is demon- 
strated in Fig. m All curves of different Js/Jb collapse 
into a single master curve, and it indicates that Eq. ^ 
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FIG. 5: Determination of the surface transition temperature 
Ts for Js/Jb = 2.0. The dashed hne is a power-law fit to the 
curve of T = 4.955. The lattice size is L = 80. 



does describe the crossover behavior during the dynamic 
relaxation. The scaling plot in Fig.[3]yields the exponents 
(j) — 0.52 and /3l^ — 0.220, as well as the special transi- 
tion point = 1.50. The crossover exponent is very 
close to the mean- field value 0.5 [H, and /S^^ and r^^ are 
in agreement with the existing results from simulations 
in equilibrium in Table |TT] and in Ref . [sll . Although the 
precision of and critical exponents obtained here are 
not very high, it is still theoretically interesting. The dy- 
namic crossover scaling form in Eq. ([5]) should be general, 
and hold in various statistical systems. 

To carry out the simulation at the surface transition, 
we fix Js/Jb at 2.0, well above r^p. At the surface tran- 
sition, where the critical fluctuation is essentially two 
dimensional, Vg and Zs in Eq. Q become V2d and Z2d- 
Around the transition temperature, the surface mag- 
netization obeys a dynamic scaling form < mi{t) >~ 
f'l3i/^,z,p(^^i/usz,^^ ^ To determine the surface tran- 
sition temperature Tg, one may search for a best-fitting 
power-law curve to the surface magnetization. Then the 
corresponding temperature is identified as the transition 
temperature Tg. We perform the simulations with three 
temperatures around the transition temperature Ts , and 
measure the surface magnetization. The results are dis- 
played in Fig. [51 Interpolating the surface magnetization 
to other temperatures around these three temperatures, 



one finds the best power-law behavior of the surface mag- 
netization at Tg — 4.955. The corresponding slope of the 
curve gives l3i/zs = 0.0570(10) at Ts = 4.955, and it is 
in agreement with the value in the 2d Ising model 
Therefore we take Ts — 4.955 as the surface transition 
temperature, which is consistent with Ts = 4.9575(75) 
obtained with Monte Carlo simulations in equilibrium 

The time-dependent cumulant U and susceptibility xii 
at the surface transition are measured, and displayed 
in Fig. [SI The slope of cumulant is 0.916(15), in a 
good agreement with 2/z2rf = 0.926(9) of the 2d Ising 
model [ig. Consistence is also observed for the suscep- 
tibility where the slope is 0.824(10), in comparison with 
l2d/ Z2d — 0.810(8) in the 2c? Ising model. We thus con- 
firm that the surface transition belongs to the universal- 
ity class of the 2d Ising model. Meanwhile, Ts — 4.955 is 
a good estimate of the surface transition temperature. 



IV. SHORT-TIME DYNAMICS ON A 
NON-PERFECT SURFACE 

In this section we investigate the nonequilibrium criti- 
cal dynamics on a non-perfect surface, i.e. Ji ^ Jg- The 
static and dynamic properties of a non-perfect surface are 
important and interesting, because real surfaces are often 
rough, due to the impurity or limitation of experimental 
conditions Q. Furthermore, the advance in nano-science 
allows experimentalists to create other structures on top 
of films artificially. We study the line defect on a surface, 
and the procedure can be generalized to other extended 
defects. 

We first consider the dynamic behavior of m2 at the 
ordinary transition. For convenience, we fix Js/Jb = 1-0. 
The profiles of m2{t) with Ji = 0.5Js, Ji = l.OJ^ and 
Ji = 1.5Js are depicted in Fig. [7] All lines look paral- 
lel to each other, and it indicates that they may belong 
to a same universality class. By fitting these curves to 
the power law in Eq. (O, we estimate P^'''^ = 0.792(18), 
0.786(6) and 0.797(33) with J, = 1.5Js, J/ = l.OJ^ and 



Surface transition 
at perfect surface 



U„(t) 



J7J|, = 2.0 T= 4.955 



Ordinary transition 
at non-perfect surface 
J./Ju= 1.0 




FIG. 6: Dynamic relaxation of the surface susceptibility and 
cumulant at the surface transition Js/Jb = 2.0 plotted on 
double-log scale. The lattice size is L = 80 and Ts = 4.955. 
The dashed lines are power-law fits to the curves. 



FIG. 7: Dynamic relaxation of the line magnetization at the 
ordinary transition on a non-perfect surface with various Ji 
is plotted on a double-log scale. The slope of the curves is 
independent of J;. 
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FIG. 8: Dynamic relaxation of the line magnetization at the 
special transition Js/ Jb = Tsp = 1.5004 on a non-perfect sur- 
face with various J; is plotted on a double-log scale. The open 
circles are a fit to Eq. (|10|l with a correction to scaling. The 
inset displays the line magnetization at J; = 1.6Ja but with 
a longer simulation time. The lattice size is L = 128 and 
T = Tc. The slope of the curves is dependent on Ji even after 
taking the correction to scaling into account. 



J; = 0.5 Js respectively. These values are consistent with 
each other and with /3°'''^ on the perfect surface reported 
in the previous section. It confirms that the defect in 
the ordinary transition is irrelevant, in term of the renor- 
malizationgroup argument. This conclusion echoes that 
in Ref. |32l |. where the impact of random bonds on 
the surface is investigated in equilibrium. According to 
the generalized Harris criterion 3j|| , defects with random 
bonds or diluted bonds on a surface are irrelevant. The 
short-time dynamic approach shows its merits in identi- 
fyin g th e universal behavior of the surface magnetization 
[s^. Tiol. m, 111]. Here we note that the line magneti- 
zation is one-dimensional, and therefore somewhat more 
fluctuating than the surface magnetization. 

Now we turn to the special transition. We perform 
simulations with various J/ at the special transition, and 
the line magnetization is presented in Fig. [H From the 
slopes of the curves, one measures the exponent (32^ / VsZs, 
and then calculates /3f = 0.260(4), 0.230(3), 0.219(5), 
0.204(6) and 0.162(3) for Ji = OAJs, Ji = O.SJ^, J; = 
l.OJs, Ji = 1.2Jfi and J/ = 1.6Js respectively, with i^., and 
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FIG. 9: Dynamic relaxation of the line magnetization with 
various Ji at the surface transition is plotted on a double-log 
scale. The lattice size is L = 80 and T — Ta. The slope of the 
curves is dependent on J; , and given in Table IIIII 



FIG. 10: Dynamic relaxation of the line susceptibility with 
various J; at the surface transition plotted on a double-log 
scale. The dashed lines are power-law fits. The slope of the 
curves is dependent on Ji, and given in Table Hill 

Zs taken as input from Table HI Obviously /Jj'' changes 
continuously with J/. 

Since there exists certain deviation from a power law 
in shorter times for the curves with a larger ratio Ji/ Js 
in Fig. [21 one may wonder whether the small variation 
in /32^ may stem from the correction to scaling induced 
by the defect line. Therefore, a careful analysis of the 
correction to scaling is necessary in this case. Assuming 
a power-law correction to scaling, m2{t) should evolve 
according to 

m2(t) = at-'3r/'^^^=(l - bt-''). (10) 

As shown in Fig. [51 such an ansatz fits the numeri- 
cal data very well, and yields = 0.258(1), 0.235(7), 
0.228(3), 0.214(6) and 0.171(3) for J, = 0.4 J^, J/ = 0.8 J^, 
Ji = l.OJs, Ji = 1.2Js and Ji = 1.6Js respectively. For 
J; = 1.6Js, we extend our simulations up to a maximum 
time t = lOOOOMCS to gain more confidence on our re- 
sults. Still 132^ varies continuously with J;, and the strong 
universality is violated. This is different from the case on 
a random surface, where the generalized Harris criterion 
states that the enhancement of the short-range random- 
ness on the surface is irrelevant at the surface transition 
in the 3c? Ising model (ssj. Our result is, however, not 
surprising, for the defect line is not a short-range ran- 
domness [31, but an extended one. As the short-range 
random surface is close to being relevant it is not 
surprising that the defect line modifies the surface uni- 
versality class. This can also be understood. The reduc- 
tion of the coupling in the defect line is somewhat like 
turning the local surface from the special transition to 
the ordinary one, and therefore gives rise to a large value 
of the critical exponent /?2^. 

To investigate the impact of the line defect at the sur- 
face transition, we fix Js/ Jb — 2.0. We measure the 
time evolution of the line magnetization at its transition 
temperature Ts = 4.955 with J/ = 0.5Js, J/ — l.OJg, 
Ji = 1.5Js and J; = 2.0Js. In Fig. [9l one observes 
that after a microscopic time tmic ^ lOOMCS, the power- 
law behavior emerges. However, the exponent (32 is Ji- 
dependent, and the strong universality is violated. This 
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TABLE III: Comparison between the numerical simulations of surface transition with a non-perfect surface and the theory of 
the two-dimensional Ising model with a defect line. Vs = !^2d ~ 1 and Zs — Z2d = 2.16(2) have been taken as input p^ . 





line magnetization 


susceptibility 


cumulant 


exponent 


P2/ZS 


(1 - 2f32)/zs 


l/Zs 




Simulation Theory 


Simulation Theory 


Simulation Theory 


Ji = 0.5Js 

Jl = l.OJs 

Ji = 1.5J, 
Jl = 2.0J, 


0.0923(36) 0.0936(9) 
0.0570(10) 0.0579(5) 
0.0301(24) 0.0307(3) 
0.0149(12) 0.0145(1) 


0.282(4) 0.276(3) 
0.356(5) 0.347(3) 
0.405(4) 0.402(4) 
0.428(9) 0.434(4) 


0.462(2) 0.463(4) 
0.475(2) 0.463(4) 
0.468(2) 0.463(4) 
0.459(9) 0.463(4) 



is similar to the case in Ref. [H, where a non- 
universal behavior of the edge and corner magnetization 
has been found at the surface transition. 

Since the surface transition is essentiaUy two- 
dimensional, one may relate this non-perfect surface to 
the 2d Ising model with a defect line without the pres- 
ence of bulk. The violation of the strong universality of 
the 2d Ising model with a line or a ladder defect is rig- 
orously proved by Bariev [i^. For the line defect, exact 
calculations show that 

2 

/32 = ^ arctan^(K;), (11) 

with 

Ki = exp(-2(J, - J)/kBT,). (12) 

The critical exponent (32 reduces monotonically, when the 
defect coupling Ji is enhanced. We measure the exponent 
(32 and compare it with the exact values obtained from 
Eqs. pTjl and (|12p. The results are summarized in Table 
mil One finds a good agreement between simulations and 
exact results. A similar behavior of the edge magnetiza- 
tion, which can be viewed as a line defect at the surface 
transition, is also observed in Ref. [11]. Our results sup- 
port that at the surface transition, the critical exponent 
(32 will change in the presence of a small perturbation. 

Finally, the susceptibility X22{t) and cumulant U22{t) 
of the line magnetization, which are similarly defined as 
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FIG. 11: Dynamic relaxation of the cumulant with various Ji 
at the surface transition plotted on a double-log scale. The 
dashed lines are power-law fits. The slope of the curves is 
independent of Ji, and given in Table HiH 



those of the surface magnetization, are also measured. 
The results are plotted in Fig. [10] and [11] Simple scal- 
ing analysis shows that X22(0 ^ ^{d-2-2P2/'^s)/zs g^jj(j 
U22{t) t'^d--^)/^'\ The estimated exponents are also 
compiled in Table IIIIl and a good consistency with the 
theory can be spotted. 



V. CONCLUSION 

With Monte Carlo simulations, we have studied the dy- 
namic relaxation on a perfect and non-perfect surface in 
the 3d Ising model, starting from an ordered initial state. 
On the perfect surface, the dynamic behavior of the sur- 
face magnetization, susceptibility and appropriate cumu- 
lant is carefully analyzed at the ordinary, special and 
surface transition. The universal dynamic scaling behav- 
ior is revealed, and the static exponent Pi of the surface 
magnetization, the static exponent 711 of the surface sus- 
ceptibility and the dynamic exponent Zs are estimated. 
All the results for /3i are compiled in Table (III Since the 
exponents Vs and Zs can be identified as those at bulk, 
it is convenient to study different phase transitions from 
the non-equilibrium dynamic relaxation. Especially, the 
dynamic crossover scaling form in Eq. ^ is interesting. 
Because of the existence of new scaling variable Jg/Jb, 
the nonequilibrium relaxation of magnetization at the 
critical temperature may not obey a power law, which 
is quite different from the general systems investigated 
so far where a power law behavior was always expected. 
This unusual nonequilibrium behavior is a consequence 
of the presence of geometric surface. 

On the non-perfect surface, i.e., with a defect line in 
the surface, the universality class of the ordinary tran- 
sition remains the same as that at the perfect surface. 
On the other hand, for the special and surface transi- 
tions, the critical exponent (32 of the line magnetization 
varies with the coupling J; strength of the defect line. 
The susceptibility and appropriate cumulant of the line 
magnetization also exhibit the dynamic scaling behavior 
and yield the static exponent 722 and the dynamic expo- 
nent Zg. The short-time dynamic approach is efficient in 
understanding the surface critical phenomena. 
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